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Phase separation of two-component untrapped and trapped Bose gases at finite temperature is studied. 
Conditions for the stability of a homogeneous binary Bose gas mixture at finite temperature and the density 
profile of the trapped binary mixture are obtained. Our extended results include the known results for zero 


temperature. 


PACS number(s): 03.75.Fi, 05.30.Jp, 64.75.+g 


I. INTRODUCTION 


Since the realization of Bose-Einstein condensation 
(BEC) in dilute atomic gases [1], a number of interesting 
experiments have been conducted to investigate multicompo- 
nent Bose gases, in which two or more internal states of 
condensates exist together in a magnetic or optical trap 
[2-6]. The dynamical response of the condensates to a dis- 
continuous change in an internal state and the dynamics of 
component separation have been studied [3]. The relative 
phase between two condensate states has been measured us- 
ing a time-domain separated-oscillatory-field condensate in- 
terferometer [4]. Spin domains in spinor Bose-Einstein con- 
densates have also been studied [7,8]. Very recently, vortex 
states have been obtained in a two-component Bose gas [9]. 

The progress in the experiments exploring dilute mixtures 
of quantum gases has stimulated intensive research on the 
properties of mixed Bose gases. There are many theoretical 
studies on the properties of trapped two-component Bose 
gases at zero temperature [10—18]. The spatial structures of 
condensates have been studied by means of the Thomas- 
Fermi approximation [10], or by a direct numerical treatment 
of the two-component Gross-Pitaevskii equation [11-13] 
and Monte Carlo simulation [14]. Collective excitations of 
two-species condensates have been studied for both the ho- 
mogeneous and inhomogeneous cases [15,17,18]. Due to in- 
terspecies interactions, the condensate mixture displays 
novel behavior not found in a pure condensate. Under certain 
conditions the system becomes unstable with respect to 
stratification, i.e., spatial separation between the two conden- 
sates, which is called phase segregation. Much of the inter- 
esting behavior of binary mixture Bose gases happens in this 
regime. For example, in a cylindrically symmetric trap, Chui 
et al. [14] and Ohberg and Stenholm [13] found the absence 
of cylindrical symmetry in the density difference of a binary 
mixture of alkali atoms in the phase segregated regime. 

In most of the calculations, the effects of finite tempera- 
ture have not been discussed. Although temperature fluctua- 
tion is included in the Monte Carlo simulation [14], effects 
of the thermal cloud on the condensate have not been con- 
sidered. How these structures change with temperature is 
thus a very interesting problem. As we know, at a certain 
temperature a Bose gas mixture may undergo Bose-Einstein 
condensation. If the two components do not interact, the con- 
densates will form and disappear independently, showing 
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two different transition temperatures. With intercomponent 
interaction turned on, the presence of one condensate will 
affect the formation of the other. This will shift the transition 
temperature, depending on the strength of the interaction be- 
tween the two components. A more important effect of the 
interaction is that it may lead to phase segregation, from 
which rich spatial structures come. 

In this paper, we pay more attention to finite temperature. 
We first derive conditions for phase separation at different 
temperatures for the untrapped homogeneous binary mixture. 
Our results agree with those obtained by analyzing the sta- 
bility of collective excitations at zero temperature. We then 
investigate properties of the trapped binary mixture of Bose 
gases by means of the local density approximation. We ob- 
tain the density profile of the condensate at finite tempera- 
ture, which includes the results given in Ref. [10] for zero 
temperature. 

The paper is organized as follows. In Sec. II we study the 
phase separation of as untrapped homogeneous binary mix- 
ture of Bose gases. In Sec. III we extend the discussion to the 
case of a trapped binary mixture of Bose gases. Section IV is 
a brief summary. 


Il. HOMOGENEOUS BINARY MIXTURE OF BOSE GASES 


Two-component ideal Bose gases are rather simple. There 
are no interactions between the components, so the two com- 
ponents behave independently. There are two critical tem- 
peratures T., and T..,, which describe the two BEC transi- 
tions. However, when the interatomic interaction is turned 
on, the two critical temperatures and the thermodynamic be- 
havior of the mixed system may change [19-22]. A new 
phenomenon, phase separation, may appear. We start with 
homogeneous Bose gases and then consider the inhomoge- 
neous case. 

The Helmholtz free energy of a homogenous binary mix- 
ture of Bose gases can be written as 


F mix =F t+ F+F! (1) 


mix 
where F, and F, are the Helmholtz free energies of compo- 
nents 1 and 2, and F” is the contribution of the intercompo- 
nent interaction. According to the thermodynamical theory 
of mixtures, the homogeneous phase is stable only when the 


symmetric matrix u, whose elements u;j are the derivatives 
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(ðuilðp;)r,y Of the chemical potentials, is non-negatively 
definite, i.e., when the determinants of all principal minors of 
ù are non-negative [23]. For a binary mixed system, the 
stability conditions are 


M1120, M0, (2) 
and 
Mit Miz 
ti | =0. (3) 
M M22 


If these conditions are not all satisfied, the homogeneous 
phase is no longer stable. The system may evolve to a non- 
homogeneous phase where one region is rich in one compo- 
nent and the other region rich in the other component. This 
phenomenon is called phase separation. As we know, F; 
takes different formulas when the temperature is below and 
above the BEC transition temperature. So, dependent on 
temperature, the behavior of phase separation will be differ- 
ent. For specificity of discussion, we assume T,.;>T... We 
now consider three different temperature regimes. 


A. The case of T=T,, 


When the temperature T>T,.,, both components are in 
the normal phase. The Helmholtz free energies can be writ- 
ten as [24] 


ut? 


y 
BF;=N, Inz; = 38s) + 2aipiN iN? i=1,2, (4) 


t 


where p;,z;=exp(6y;), and \;= \2 ah? /m;kT are, respec- 


tively, the density, fugacity, and thermal wavelength for the 
component i, and a; is the s-wave scattering length. The 
particle number N; is given by 


V 
Ni= 3 832i). (5) 


L 


The expression for the contribution of the interaction be- 
tween the components, which is valid for a homogenous 
state in all temperature regimes, takes the form [25] 


BF! =a (A7 +AS)N N/V, (6) 


where a; is the scattering length between the two compo- 
nents. From Eqs. (4) and (6), we obtain the chemical poten- 
tials as 


Bui=BuP+4ai pii tanp +43), (7) 


Ber= Bey? + 4arp 3 tanpi (AF +AS), (8) 


where po are the component chemical potentials for the 
ideal Bose mixture when the interaction is turned off. The 
second terms come from the mean field of atoms of the same 
kind and the third terms from that of atoms of the different 
kind. 


PHYSICAL REVIEW A 61 063613 


We now analyze the stability of the binary mixture. For 
convenience later on we use By instead of 2, but still keep 
the same notation 4. From Eqs. (7) and (8), the diagonal 
elements of & are 


ABE, OBES” 


+4a,d? 
Ki api ai a,Qj, (9) 
dB, dbus? A 
= = +4a5N35, (10) 
Ha ap dBpr : 


and the off-diagonal elements of 2 can be written as 
H= My =(Ai+A5) ap. (11) 


By taking derivatives on both sides of Eq. (5) with respect to 
pi, the derivative of the chemical potential pO for the ideal 
Bose gas is found to be 


Be dj 


Pi gueza AFP) 


where we have used the relation zdg,(z)/dz=g,_,(z). The 


(12) 


determinant of the symmetric matrix j can then be written 
as 


dj 
+4a;r? 


A=det(z)= I] - 
i=12 | galez (3 p;)) 


—(AG+N3)7atp. (13) 


Let us consider the two limiting cases of T— and T 
—T.,. When T—~™, then \—>0, which leads to 


8 2832 AFPA P; . 


We then have w,,;=1/p,, M22= 1/p2, and A ~=1/p;p2. These 
are all positive. Thus, at high temperature the homogeneous 
phase of the binary mixture is always stable, as it should be. 
No phase separation happens at a very high temperature. 
When T—T.,, from Eq. (5) \}py—g3p(1)~2.612; 
hence 
81283207 P= 8 nL). 


In this case, 41; and uz reduce to 


Mi1=4a1d2,, (14) 


A2 
27 = 
& iol Sap (2.612(T.» ITa y) 


+4ayd5, (15) 


and 

M2 
‘eee CINTAT N 
Ac HAS) ai (16) 


Ac =4a A? +4ayr5 
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where \,; denotes N; at temperature T=7,;. When the in- 
teratomic interaction of the component | is attractive, 411 
<0. Thus the homogeneous phase is never stable for attrac- 
tive Bose gases, as is well known. In the following discus- 
sion, we shall consider only the case when the inter- and 
intracomponent interactions are repulsive. From Eq. (3) we 
obtain the stability condition for the homogeneous phase, 


^o /a3 ' 2 Ge (17) 
Elg za (2.612(T o/T a) az’ 


where 


m +m 
o my 2 (18) 


all 2Vm\m, 


A threshold value aj, of a; is determined by A,,=0 to be 


1/2 
à2/a3 Vaja 


=] wa 14 
8 1283/9 B31 Teal T a) X 


u 
a2 


(19) 


When a j)>a{, phase separation will occur at a temperature 
above T... 

As a simple example, when 7,.,=T..=T,, the condition 
(17) reduces to 


Ax4a,a,— x7ai,>0. (20) 


Since g,(z) is monotonic, A will increase with temperature. 
So if condition (20) is satisfied no phase separation occurs in 
this system at T>T.,. If 4a;a,< X dio, the homogeneous 
phase will not always be stable at T>T,.. There is a tem- 
perature below which the homogeneous phase becomes un- 
stable and phase separation takes place. For the simple case 
of a}=a3, m;=m3, and p; =p, from Eq. (13) the critical 
temperature for the phase separation is determined by 


Ni 


=2a;2 4a). (21) 


8108 3201p 1)) 
Generally, the critical temperature for phase separation has 


to be found numerically. 


B. The case of T,..<T<T, 


At a temperature between the two critical temperatures for 
BEC, component 1 has undergone condensation, while com- 
ponent 2 is still in a normal state. The Helmholtz free energy 
for component 1 should be written as [24] 


V Z472 2n72 
pry = eee eae dik 


i 
T\32 | 7 \3 
= = | 


(22) 


aN 


2y72 
H1 
V 


1+2 


4 (1)+ 
=~ 3 85/2 
Mi 
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where the relation N,;,=N,[1—(T7/T,,)°*7] for the number 
of atoms in the condensate has been used in the derivation. 
The Helmholtz free energy of component 2 is still given by 
expression (4). Including the contribution from the intercom- 
ponent interaction, we have the chemical potential 


T 3/2 | T 
Tei | Tei 


The elements u4; and u; of matrix ú can then be calculated 


as 
T 3/2 T 3 
142( 2 | ata ait 


IB My 
J 1( Mi +3), (25) 
P2 


3 


1+2 +apa(A?+AS). 


(23) 


Bei=2a,djp1 


_ IBM 
Mit a 


= Iai 
Pi al 


M127 


while up is the same as Eq. (10). Finally, we obtain the 
determinant A of matrix %4: 


T \32 T \3 
A=2a,\7/14+2 = 
; Ta Tea 
M% +4a vajaga; 
EE 22 1TA2) 4in2. 
gin(8g32 A2p2)) 


(26) 


When T=T,,, A reduces to the same expression as Eq. (16). 
At T=T,», determinant A becomes 


ra) ~ | ra) 
Tai Ta 


= (AHAS) ai. (27) 


Aa=8a4 aA N2] 1+2 


Another threshold value of a; for phase separation is 


ra 3/2 | ra] ] 1/2 
2 = . (28) 


da \2a,a3 
3 Tz Ta 


2 


a 


When the scattering length a, satisfies al,<a 12<4i2, phase 
separation occurs below a certain temperature in the interval 
(T.2,T.,). When the two critical temperatures coincide, i.e., 
T1 =T 2, the two thresholds become identical as 
ah=ah=2vVa;ax, 

which is consistent with Eq. (20). The critical temperature of 
the phase separation for the binary mixture with al,<a 12 
<a can be determined from A =0. 


C. The case of T<T.2 


When the temperature T<T.2, both of the two compo- 
nents have undergone Bose-Einstein condensation. The 
Helmholtz free energies for components | and 2 are of the 
same form as Eq. (22). The chemical potentials u; and u3, 
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T/a 


0 4 8 12 
S 
FIG. 1. Dependence of the critical temperature of phase separa- 
tion on s=x%a;2/ Vaaz. Here we take N, /a1=10, Nc2/a2=10, 
and Ta =2T.2. 


and the elements of matrix ú then take corresponding forms 
similar to Eqs. (23)-(25). It is easy to derive the result that 


T 3/2 T 3 
142 | [5] -ajta 


= 2 
A= [I 2a,X7 F F 
(29) 


i=1,2 


Phase-separation conditions are determined from A =0 as 
in the above cases. At temperature T=T.2, this leads to the 
threshold af, for a). When a,.< os phase separation oc- 
curs only at T<T,2. Another threshold aĵ, for a4) corre- 
sponds to the limit T—0, where 


Ay~4ayaohth3—(ATHN5)2ai0, (30) 
which implies 
aty=X Va) a. (31) 


When a{,<a igus phase separation occurs at a tempera- 
ture below some T, € (0,7). If ay2<a{,, no phase separa- 
tion occurs in this temperature regime. For mı =m, expres- 
sion (31) reduces to a°,= aa, recovering the stability 
condition at zero temperature derived in Ref. [16], based on 
the stability analysis of elementary excitations. 

For a{)<a).< af, the temperature of phase separation is 
determined by 


elz) r 


which, as a quartic equation in T°”, admits an analytic solu- 
tion. For the simplified special case of To =T 2, the result is 


2,2 
X aiz 


ajaz’ 


(32) 


reti -yanina Ta (33) 


The above discussion shows different behaviors of the 
phase separation in different temperature regimes. In Fig. 1 
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we plot the critical temperature for the phase separation as a 
function of s= ya 2/ Vaia», setting the other parameters at 
Noy /4,=No2/a,=10 and T. =2T.2. Results near the two 
BEC transition temperatures T, and T, are not very accu- 
rate according to the discussion in Ref. [24]. 


IHI. TRAPPED BINARY MIXTURE OF BOSE GASES 


In real experiments, systems are not homogeneous. Bose 
atoms are trapped in an external potential. We shall extend 
the above discussion for the homogeneous binary mixture to 
a trapped binary mixture, using the local density approxima- 
tion. 

In the discussion below, by T, and T.) we mean the 
BEC transition temperatures of the two trapped components. 
We still assume that T.;>T... When T>T,.,, both compo- 
nents are in the normal phase. As given in Ref. [26], the local 
grand potential is 


A)T È [Aj gs2) + 2aipi N] 
+an(AÏ+AS)PIP2, (34) 

and the density distributions 
pi(r)~y 83E), (35) 


p(¥)~Az 832l Er), (36) 


with 


¿=z expl- BV(r) —4a,\jp1 (1) — 2a 4( AT +23) p2(r)], 
(37) 


&=z> exp[ — BV(r) —4ap\3p9(r) — 24 9(A7+A3) pi (0). 
(38) 


The chemical potential may be determined from the normal- 
ization relations 


~i=f dr p,(r), (39) 


w= { dr p2(r). (40) 


For the general power-law potential 


p l 


tE 


x q 


Lı 


7 
A 


L; 


V(r)=€; +63 


> 


Ly 


with e;, Li, p, l, and q being parameters for the potential, 
by integrating the local grand potential over the whole vol- 
ume, we find the total grand potential of the system to be 
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O nix= > 


—2a;r;! Z; 
F745 pagr als nC) 24d; H3232,9(Zi Zi) ] 


Wa yo(j +03) 
8310312, 921 222): (41) 
Leer 7 
where 
1 1 1 1 a xyi 
Q= tS tt Ay )= > 


ij=l hd CoD ea 


_ 8L,L5L31(p.l,q) 
eyes et” 


? 


I(p,l,q)=(plq) ‘PO /p)POM)P (1/4). 


Using a procedure similar to that of Refs. [20,21], we can 
discuss the thermodynamic properties of the binary mixture, 
and then examine the stability of the thermodynamical states. 

A more interesting case is when 7..<T<T,,. Since com- 
ponent 2 is in a thermal state, its density distribution can still 
be written as 


palt) ~ )3 ° gzl Er). (42) 


However, the density distribution of component 1 has no 
simple unique form over the volume. In the domain without 
a condensate, the density distribution can be characterized by 
the parameter mı, which is similar to the case when T 
>T,,. In the domain with a condensate, the density distribu- 
tion is quite different. From Eqs. (6) and (22), in the conden- 
sate domain the Helmholtz free energy density can be written 
as 


BIi=-A  8sa(1)+2a M pia pit BVE)p: 
halt pipa, (43) 


where p}, is the condensate density, which, together with the 
uniform thermal density distribution pio, gives pı(r) 
=P;ı;(r)+ pọ. The chemical potential of component 1 then 
takes the form 


ui=G(pi(r)+ pio) + G 2p2(r), (44) 
where 
4ma;h? 21a jyh?(m, +m) 
Gi= > Gy= ale ; . (45) 
m; mm 


Thus, we have 


pi(r)=[ui — V(r)—Gip2(r)])/G11— P10- (46) 


The density distributions can be characterized by the two 
chemical potential parameters u; and u3, which are fixed by 
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relations (39) and (40). Once the chemical potentials are de- 
termined, all other thermodynamic properties can be calcu- 
lated. 

As we already know, the critical temperature for a single- 
component Bose gas in the general power-law potential is 
T.~(N/m??) "TD [20,22]. For the special case of N, 
<N; in the binary mixture, Te2&T.1 is satisfied. When the 
temperature T is close to T..,, almost all the atoms of com- 
ponent | are in the condensate. Neglecting the contribution 
from the thermal atoms, we have the density profile of com- 
ponent | as 


pir)~= [ui = V(r) VGi1, (47) 


where u; is fixed by expression (39). Atoms of component 2 
are then trapped in an effective potential of the bare trap and 
the interaction field of component 1: 


[ui = V(r) ]Gio 
Gii 


Vig7= Vn) 4 for r<ro, (48) 
where rọ is the border of the condensate of component 1. 
When G}.>G,,, the minimum of this effective potential for 
component 2 is no longer at r=0, but at r=rg. As we know, 
condensation first takes place at the minimum of the poten- 
tial. The condensation of component 2 may thus first appear 
at the surface surrounding the condensate of component 1. 

When T<T.», every component has two domains: the 
condensate domain and the thermal domain. The condensate 
domain consists of the saturated gaseous part with density 
Pa- eel) and the condensate part with density p;, 
=P; Pio. AS a consequence, there may exist three kinds of 
subdomains: (a) with both components 1 and 2 in the thermal 
states, (b) with one component in the thermal state and the 
other in the condensate, and (c) with both components in the 
condensates. In subdomain (a), the density distributions are 
similar to those for T>T,,. Subdomain (b) is like the mix- 
ture at T.9<7T<T,,. In subdomain (c), the free energy den- 
sity can be written as 


Bfi= -N Pgsp(1) + 2a;Fp;— aj; Pist BV() 0}, 


(49) 
and the interaction part as 
Bf'=a (A+A) Pipa. (50) 
The chemical potential is then given by 
M1 =G (pi (4) + pio) + Gi2p2(r), (51) 
b= Ga9(p2(r) + pao) + G 1291 (4). (52) 


The density profiles of thermal atoms in the condensate do- 
main are position independent, so the variation of the total 
densities comes from the variation of density distribution of 
the condensate. We rewrite Eqs. (51) and (52) in the forms 


M1 =Gy1P1s(0) + Gy2p2,(4), (53) 


[y= Gaopos(¥) + G i2P1s(£), (54) 
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where 


M1 = M1 — 2G pio G 12P20> f= a= 2G P29 G 12P10- 


From Eqs. (53) and (54) the condensate density distributions 
may be derived as 


= HiG- MG 1 + G> Gy 


Pis V(r), (55) 
GuGn-Gh GuGn- Gh 
[LG 11—- uG Gy-G 

go 117 Mi a 12 - V(r). (56) 
G\;Gy— Gi GiG- Gj, 


The above expressions for the density profile are similar to 
those obtained in Ref. [10] for zero temperature. At zero 
temperature these expressions recover exactly the earlier re- 
sults. The denominator G,,G5.— Gi is the quantity charac- 
terizing phase separation at a finite temperature. It is the 
same quantity that determines whether there is a phase sepa- 
ration in trapped binary Bose gases at zero temperature [10]. 
A trap makes condensates at finite temperature more like 
condensates at zero temperature. This is one reason why 
theories for trapped Bose gases at zero temperature can de- 
scribe experiments well, although experiments are not con- 
ducted at zero temperature. 


IV. SUMMARY 


In the above, we have investigated untrapped and trapped 
two-component Bose gases at finite temperature. Contribu- 
tions from the thermal cloud atoms have been taken into 
account. By analyzing thermodynamic stability, we have ob- 
tained conditions for the phase separation at different tem- 
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peratures for the homogeneous binary mixture that agree 
with those for zero temperature found by other authors. 
These results show different behavior of phase separation at 
different temperature regimes. We have extended the discus- 
sion to trapped binary mixtures of Bose gases by means of 
the local density approximation. Expressions for the density 
profile of condensates have been derived. It is found that 
phase separation may occur in traps even at a temperature at 
which no phase separation would occur in the homogeneous 
case without trapping. This kind of spatial separation comes 
from the effects of trapping. The trap focuses condensed at- 
oms in the center region where the condensate is close to the 
homogeneous condensate at zero temperature. 

Recently, the dynamics of phase separation have been 
studied in an experiment [3]. The time scale of phase sepa- 
ration has been estimated by studying the fastest growing 
modes [27,28] in homogeneous binary Bose gas mixtures at 
zero temperature. As shown in this paper, the behavior of 
phase separation can be modified by the existence of the 
thermal cloud. Now a similar question arises: How does the 
thermal cloud affect the time scale of phase separation? The 


determinant of the symmetric matrix a may play a similar 
role to the fastest growing modes at zero temperature. When 
the determinant becomes negative, phase separation occurs. 
The dependence of the time scale of phase separation on 
temperature may be deduced from this quantity with some 
assumptions. Of course, the time scale will also be affected 
by the presence of a trap. This is under study. 
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